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Abstract 

We focus on the long time behavior of complex networks of reaction- 
diffusion (RD) systems. We prove the existence of the global attractor 
and a I/°°-bound for a network of n RD systems with d variables each. 

This allows us to prove the identical synchronization for general class of 
networks and establish the existence of a coupling strength threshold value 
that ensures such a synchronization. Then, we apply these results to some 
particular networks with different structures (i.e. different topologies) and 
perform numerical simulations. We found out theoretical and numerical 
heuristic laws for the minimal coupling strength needed for synchroniza¬ 
tion relatively to the number of nodes and the network topology, and 
discuss the link between spatial dimension and synchronization. 

Reaction-Diffusion systems, complex networks, attractor, synchronization. 


1 Introduction 

Networks of dynamical systems appear naturally in the modeling of numerous 
applications. The simplest mode for the coordinated motion between dynamical 
systems is their identical synchronization when all nodes of the network acquire 
identical dynamical behavior. Such cooperative behavior has been observed in 
natural or artificial systems such as neural networks, chemical and biological 
systems, computer clocks, social networks... In this paper, we will focus on 
RD systems networks that can be seen as neural networks. Besides, a classical 
question in dynamical systems is the existence of the attractor: basically, a set 
that attracts the trajectories for large time. 

Contributions. There are three main contributions in the present paper. 
First, we prove of the existence of the network attractor, and therefore, within 
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this attractor, we analyze the synchronization behavior. Finally, we perform 
numerical simulations with different kind of initial conditions and study numer¬ 
ically the spatial effects on synchronization and pattern formation. 

Mathematical framework and preliminaries. First of all, we introduce 
the mathematical framework we will use throughout this paper. Mathematically 
speaking, the network is represented by a graph, the nodes of which are a d- 
dimensional RD system and the edges correspond to the coupling functions 
between these subsystems. The general system reads as: 

Uu=QAU,+F(U,} + H,(Ui,...,Un),iG{l,...,n}. (1) 

In this equation, each variable Ui represents a function from fl x IR+ into 
n is a bounded domain of and F : —>■ is the nonlinear reaction term. 

For all i G {l,...,n}, Jli : —)■ is the coupling function between nodes 

whereas Q is a diagonal matrix of with positive coefficients. If we add 

boundary conditions to Q, we obtain a general reaction-diffusion system. We 
will not go into details concerning the existence of the semi-group of Q. We 
refer to [191 |20l |25l |29] or [151 [171 |22j |^ |^ , for classical results on the existence 
of semi-group in or in spaces. Our first theoretical result is the 

proof of the existence of the global attractor for a particular class of networks of 
type 0 that generalizes the FitzHugh-Nagumo (FHN) equations. Recall that 
(FHN) equations are a simplification in two variables of the HodgKin-Huxley 
model of four equations for the action propagation in nerve, m [IH1I23]. A 
good qualitative analysis of the (FHN) reaction-diffusion system is given in [26] . 
while in [5] we gave a first analysis of a particular network of (FHN) reaction- 
diffusion systems. Here we extend some results of [5], and also of [H] where 
a system of two scalar equations, in which the diffusive term appears only in 
the first equation but not in the second one, was considered. We present some 
results for a network of n partially diffusive systems with d equations. Indeed, 
we suppose that we can split the system 0 into two subsystems, diffusive and 
non-diffusive, with s and d — s equations. Therefore, we set for all i € {1,..., n}, 
Ci = {ui,Vi), and write Q in the following way: 

f Uit = F{ui,Vi)QAuiHi{ui, ...,Un), on Hx]0,-|-oo[, ze{l,...,n} 

( Vit = —tT(a;)ui-I-<I>(a:, Ui) on Hx]0,-|-oo[, 

^ ( 2 ) 

with Neumann Boundary conditions on dfl, and where Ui take values in K®, 
1 < s < d whereas Vi take values in Q is a diagonal matrix in M®^® 

with coefficients q\j G {1,..., s}, and take values in K®. We use the classical 
notation ut for ^. This means that diffusion and coupling terms appear only in 
the s first variables of each subsystem of the network. Finally, a{x) is a matrix 
in u(d-s)x(d-s) verifies: T,fZi crjiix)viv^ > ^^^=1 for all k G 

N, k odd, and some positive constant a, and with bounded derivatives. The 
application $ takes values in Under some conditions on functions the 

system ([^ generates a semi-group on Jif = (L^(H))”'^, see [HH^j^Sj. Before 
going into details of the analysis of system ([^ , we present some key features for 
systems with one variable and two variables used to prove the existence of the 
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global attractor. These techniques will be generalized to system ([^ in section 
Let us start with the following equation: 

Ut = Au — + u, (3) 


considered in a bounded domain Q with Neumann boundary conditions. Mul¬ 
tiplying (§ by u gives, 

^ + 2 [ |Vup = - [ u^) (4) 

dt JQ JQ JQ 

<- 6 [ + K, (5) 

Jn 


for some constants S and K. By Gronwall Lemma, there exists a constant K' 
such that. 


<K', 


for all initial conditions in L^{fV) and for t large enough. Now, integrating ([^ 
between t and t r for a given constant r gives. 



|Vitp < K for another K. 


( 6 ) 


Multiplying ([^ by by analog computations, we find that there exists 

a constant K" such that 

/ < K", (7) 


for all initial conditions in and for t large enough. 

Also, multiplying ^ by —Au gives. 


d 

dt 



= — 2 f (Au)^ ^ — u^)Au 

Jfi Jq 

< — / (Au)^ + h / u*) by using Young inequality. 

Jn 2 Jq 


Therefore, thanks to ([^ and Q, we deduce, by using uniform Gronwall Lemma 
(see appendix), that 

I |Vup<iL, 

Jn 

for a constant K and for all initial conditions in for t large enough. 

This gives the compacity of trajectories of (§ thanks to the compact injection 
of in . Now, we consider the system with two variables. 


Ut 

Vt 


Au — u^ u'^ uV 
—Sv u 


( 8 ) 
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Multiplying the first equation of Q by u and the second by v, integrating, using 
Green formula, Young inequality, and Gronwall lemma leads to: 


[ {v? + < K, 

Jn 


for a constant K, for all initial conditions in x L^(n), and time large 

enough. Then, we will show the same result in x for all k gN*. 

As we saw above, this result is true for fc = 1. We multiply first equation of (|^ 
by and the second by sum the two equations and integrate. We 

obtain: 


_d 

dt 


^ I / i,2fe-2|yy|2_ / (y2fc+2^^2fc+1^^2fc)^ / -6 I V 


2k 


Using Young inequality ab < provides: 


^ / u^+ 72 / u"''+"+73 / w"''+74 / v^'^-6 / u 

™ Jn Jn Jn Jn Jn Jn Jn 

with 72 < 1 and 74 < 6. Then, since < 2k, there exists constants 7 and 
K such that: 

d 


,2k 


^ I < _ / (y2fc^^2fe) 

dt ,1a .la 


K. 


It follows that there exists a constant K (depending on k) such that: 

[ +^2'=) < K, 

Jn 

for all initial conditions in L?{fl) x and time large enough. By analog 

computations as for (§ , we obtain: 


[ |Vu|2 < K. 
Jn 


It remains to consider |Vup. By multiplying gradient of the second equation 
of ® by Vu, we obtain: 


it follows that, 


which gives. 


7 ^ 

— / iVuP =-2(5|Vz;P+ 2Vu.Vu, 
at Jq 

I |Vu|^ < -S\\7vf + K, 

Jn 

[ \^v\^<K, 

Jn 


d 

dt 
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for t large enough. We have bounds in L'J for all q G [1, +oo[. We can obtain 
bounds in L°°{n) thanks to a result in JT7\ . 

Note that we can apply the same techniques for the generalized system: 


( ut = Au-\-f{u,v) 

\ Vt = —Sv-\-g{x,u) 

if the following conditions hold 

uf(u, v) < — 6 iu^ + S 2 UV + i53 

and, 

\f{u,v)\ < Si\u\^~^ + 52 |u| + 63 

\^{x,u)\ < c\u\ and \^{x,u)\ < K, (12) 

where p> 2 , Si > o,ze{i,...,3}. 

Paper Organization After the present introduction, the theoretical results 
are presented in section 2 and 3. In section 2, we prove, under some assump¬ 
tions, the existence of the global attractor for the network of FHN-type ([^ 
in (L^(n)‘^)”, for any network topology. This allows us to show, in section 3, 
theoretical results on the synchronization onset of system with linear cou¬ 
pling functions. Then, we apply these results to complete and ring networks. 
In section 4, we present the numerical simulations for fully connected and uni- 
directionally coupled ring network. Each node, in the graph, is represented by 
a system of FitzHugh-Nagumo reaction-diffusion equations. This gives an in¬ 
sight on the relation between the number of neurons and the minimal coupling 
strength needed to reach the synchronization, with a particular attention on 
the effects of spatial dimension. For fully connected networks, our numerical 
simulations show that the minimal strength value for synchronization follows a 
“i” law, independently of the patterns induced by initial conditions. In unidi¬ 
rectional coupled ring network, the minimal strength value for synchronization 
follows a “n^” law. Our conclusion is left to the last section. 


(9) 

( 10 ) 

( 11 ) 


2 Existence of the global attractor 

Now we prove the existence of the global attractor for the dynamical system 
§ in ^ = (L^(f2)‘^)". The global attractor is a compact invariant set for 
the flow that attracts all trajectories (see for example n mi US mi)- The 
existence of the global attractor is essential since it is a set where the solutions 
asymptotically evolve. In particular, all the patterns that we will see later in our 
numerical simulations, belong, for enough large time, to the global attractor. 
Also, our proof of synchronization given in section uses L°°-bounds that we 
prove in this section. Now, we specify some assumptions that we will assume 
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throughout the article. First, we assume that for all i € ,n}, and for all 

j € 


d—s 


< -S,\r4\P + S2\ui\J2 \<r + + <^4, (13) 

k^l 

with p > 2 , 61 , 62,63 > 0 , 0 < Pi < p — 1 , and, 

s d—s 

\F^{u.,v,)\ < s,\uir^+ 62 Y.\u^r + i^'i + ^4- (w) 


k=l 


k=l 


Condition (131 generalizes (101. It indicates a decrease of order p at infinity 
and permit to obtain bounds in L‘^ spaces. Condition (14) generalizes (11) and 
allows us to apply Young inequalities in order to obtain bounds in Ff^. In our 
case, a typical example for which (13)-(14) hold, is given by a function F where 
the component j reads as: 

p—2 s 


F^ul...,ut,vl,...,vf '^) = -ap-i{ul)P ^+Y^ E n(^^)“'‘ 

k—Oaki + ...+ocks—k 1 — 1 

d—s 

1=1 

with ttp-i > 0, G ffi and p even. This simply means that, relatively 

to Ui, F^ is polynomial of several variables, with the dominant term given by 
with negative coefficient, and p even. The other terms have a degree 
lower than p — 1. Whereas F^ is a linear function of Vi. 

Moreover, in order to maintain the effect of the decrease condition (13), we 
suppose that the coupling functions have a polynomial increase lower than p— 1. 
This reads as: 


\Hl{ui, ...,Un)\ < <^ 4(1 + E 0 <pi < p - 1 . 


(15) 


k=l 


Finally, we suppose that for all j S {1,..., d — s}. 


— [x,Ui) 

OXk 




(16) 


and. 




< ^5. 


(17) 


Conditions (16) and 0 generalize condition ( |12[ ) . They are not very restrictive 
and include functions $ with spatial heterogeneity, that allow rich behavior. 
















Complex network of reaction-diffusion equations 


7 


bifurcations and pattern formation (see mil). We deduce from ( [16| and dl 
that for all j G {1 ,d — s}, 


\^^{x,Ui)\ < ^ 6(1 + ^ ( 18 ) 

Let us remark that all these assumptions appear naturally in the proof of the 
existence of the attractor of system . They appear also in mm- 
The following theorem gives the existence of the global attractor. 

Theorem 1. Under assumptions ([^-([^, the semi-group associated with ([^ 
possesses a connected global attractor in ^ Furthermore, 

is bounded in 

The proof of theorem m relies on a general result that gives the existence 
of the global attractor in Banach Spaces, see 153]. If there exists a bounded 
absorbing set tfd in which means that verifies the following condition; 

for all bounded set B C J^,3tB',yt > tB,S{t)B C (19) 


and if, 


for all bounded set B C Jtf,3tB;Ct>tBS{t)B is relatively compact in 

( 20 ) 

then the w-limit set of tfS, is an invariant connected compact set that attracts 
all the trajectories. Therefore, we divide the proof of theorem [T] into four parts, 
which for the reader’s convenience, we present as different lemmas. Before 
going into details let us briefly present the sketch of the proof. We first show. 


in lemma m the existence of a bounded absorbing set in Jif, that is (19). Then, 
in lemmaM we prove a result of compacity for trajectories, that is ( |20[ ). More 
precisely, we establish the existence of a bounded absorbing set in (fl^(n))"'^. 
The result follows from the compact injection of in LF‘{fl). Note that 

proving lemmas m and gives the global attractor existence. After, in lemma 
I) we obtain the (L“(fl))"‘^-bound: we show (L9(fl))”'^-bounds for all q G N. 
Finally, we prove theorem m by using a result that links and L'J-norms for 
linear parabolic equations, see m- 
Let us introduce the following notations: 


l^lp.O = 


L\^dx) ^ , 


( 21 ) 


if u is a real or vector valued function. We also use 


to denote the euclidian norm for a real vector. The lemma below establishes 
the existence of a bounded absorbing set in Jif. 
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Lemma 1. There exists an absorbing bounded set in Jff, that is, there is a 
constant K, such that for all initial conditions in : 

{W\l,n + \v\ln){t) < K fort large enough. 


Proof. The proof mainly relies on the presence of —5i\ul\‘^ in (13) and the 
condition on cr(a;). In the following, parameters S, c and K are generic constants. 
We multiply each scalar equation in (i), for all i and j, by uj or vl, and sum. 
We show, by using Young inequality, a6 < and as it has been done 

for equation ([^ that: 

+ I^l2,n) + ^i(l“l2,n + kl2.n) + + ^ 3 |Vu| 2 _n < S 4 . (22) 

Then by Gronwall inequality we have for t large enough: 

(kiln + klln)(k < (23) 

□ 


Now, we will show the compactness of trajectories in Jif, by establishing 


(20). We have: 


Lemma 2. There exists an absorbing bounded set in (iJi(n))"'^, that is, there 
is a constant K, such that for all initial conditions in and t large enough. 


|Vu(t)|2,n < K. 


Proof. By integrating (22) between t and t + r, we have for t large enough and 
Vr > 0: 


nt-\-r 


|Vu 


2 ,n + ^ 


nt-\-r 


kip,a <K-k 6 r. 


(24) 


Now, we multiply each component of the first equation of (§ by —Aul, we 
integrate and sum over i and j. We obtain: 


1 d 

2 dt 


|Vw|ln = -^^f / {F^u„Vi)Aul + qjAulAul + ...,Un)Au{\ , 

i=ij=i / 


(25) 


which, thanks to (14) and (15) leads to: 


- J IL {3 n O LL — * IL 

2^|vkln+g|Akia<cEE / (i + E+ E+ EKr)I 

. 1=1 k=l 1=1 


i=l3=1 • 
n s 


s'YY / (5^(1 + YI“!r'+ 1 : 14 i +1:1 

i=l 3 = 1 ^ 1=1 k=l 1=1 


JIPI) 
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where q = qi and c is a generic constant. It follows that: 

n s /. s d—s n 


2-1J-1' 


/n 


z-i 


k^l 


Z -1 


< c(l + \u\Ip_I q + 1^12,n)- 


Thanks to the techniques we use in lemma we can prove that for t large 
enough: 

\u\Zzln < K. (26) 

Then we can apply the uniform Gronwall lemma (see appendix), and show that: 
|Vu|f f 2 (f) < K for t large enough. 

It remains to find a bound for |Vu| 2 ,n- For all i S {I,..., n}, and k S {I,..., N}, 
we have: 


d 

2dt 


Jn 


(27) 

with Vix^. = ■ Using Young and Cauchy-Schwarz inequalities, we find, thanks 

to (16) and (17), for enough large time: 


^kia:fc|2,n < ~ 2 


(28) 


Finally, the result follows by using Gronwall inequality and summing over i. □ 

Lemma 3. For all g € N, there exists an absorbing bounded set in (L'^(r2))”‘^, 
that is, for all q G N there is a constant Kq such that for all initial conditions 
in Jif : 

(.Htn + M t large enough. 

Proof. In the following, TT is a generic constant. We multiply the first equation 
of © by ‘^ul, for all i,j and integrate. We obtain: 


1 d 
2k dt 


4^= / 

Jq Jq. 

(29) 


But we have: 


/ q^Z^ul\u’f^ ‘^ul=-{2k-l) I \Vul\4u: 
Jq Jq 


t|2LJ|2/c-2 


(30) 


Thanks to (13) and to (15), it follows that: 

d—s 


1 d 
2k dt 


Jq Jq 


I fclpi. J|2fe-2^J 


( 31 ) 






Complex network of reaction-diffusion equations 


10 


Hence by Young inequality, 



'^~j^W\2k,n < 2fe,n + 

(32) 

By Gronwall lemma: 

VI 

q 

(M (M 

(33) 

for t large enough. 

Besides, by similar techniques, we have: 



2fc.n — “O-ki 2fe,n + 

(34) 


Then, the result follows by Gronwall Lemma. □ 

Proof of theorem^ It remains to prove the L°°-bound . For all i G {l,...,n}, 
and for all j G {1, we have: 


Jo 

where ^ represents the semi-group associated with Aip -|- (/j = 0 and 

Neumann boundary conditions. We know, see m, lemma 3 p 25, that ^ 
verifies: 

|■^(^)<^|oo.n < cm{t)~^e~^^(p\2N,n, (36) 

where m(t) — min(l,t), A is the smallest eigenvalue of the operator I — A, 
and c is a positive constant. This allows us to conclude. □ 

3 Identical synchronization of a network of n 
reaction-diffusion systems 

Now, as we have shown the existence of the global attractor, we can think about 
synchronization of solutions within this attractor. Therefore, in the present sec¬ 
tion, we focus on this ubiquitous phenomenon and restrict ourselves to identical 
synchronization, in networks with linear coupling. We will determine sufficient 
coupling strength values ensuring the identical synchronization in complex net¬ 
works 1^. We refer to 0 m uni in m im for relevant works on the synchro¬ 
nization phenomenon. These values reveal the dependence on the number of 
nodes, the coupling configuration, the properties of each subsystem. Such a 
result holds in case of spatial heterogeneity. This point will be discussed in 
more details in the last section. Technically, to establish the synchronization, 
we exhibit a Lyapunov function for the network. Since we obtained L°°-bounds 
in the previous section, it is always possible to find such a Lyapunov function 
provided that the coupling strength is large enough. We rely on previous works 
for networks of ODE’s (see [6l|7]) and assume that the connectivity matrix has 
zero row and column sums. However, we extend the results to networks of 
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reaction-diffusion systems and improve computations, which lead us to find the 
same threshold synchronization values, see below and theorems and 

Definition 1 (see [2]). Let U{t) = (Ui(t),U 2 (t), ...,Un(t)) be a network with a 
given topology. We say that U synchronizes identically if, 


lim \Ui{t) - Ui+i{t)\ 2 ,Q = 0. 

>-+oo 


2=1 


where notations are those of ([^ and (21). 
We consider the network, 


Uit = F{ui,Vi) QAuiie n} 

Vit = -a{x)vi + ^{x,Ui), 


(37) 


with Neumann Boundary conditions on dfl which is a particular case of (§ with 
Ui = {ui,Vi), where : 

n 

Lfi{u\, -, Uji) = ^ ) CikUk- 

fc=l 

We will assume that the matrix G = {cik)i<i.k<n has vanishing row and column- 
sums and non-negative off-diagonal elements, i.e., Cik > 0 for i ^ k, and 


Cii — ^ ^ Cik — ^ ) Cki- 

k—l,k^i k—l,ki^i 


(38) 


The connectivity matrix G defines the graph topology as well as the coupling 
strength between nodes. Indeed, there is an edge between node i and node k if 
and only if Cik > 0. We consider here an arbitrary connected and directed net¬ 


work of linearly coupled RD systems satisfying (38). Obviously, symmetrically 


networks with vanishing row-sums are a particular case of the one considered 
here. In order to ensure the existence of the global attractor, we assume that 
the network (37) verifies the assumptions (1T31)-(1T7|- Note that ( |I^ holds au¬ 
tomatically thanks to the linear coupling. Before going more into details, let 
us summarize the main ideas of the section. We establish the synchronization 
result in theorem which relies on the statement of a Lyapunov function V that 
reads as the sum of the norms of all vectors Ui — Uj. More precisely, the first 
step exhibits a diagonal definite positive matrix A that annihilates the effects 
of nonlinear terms of (37), this is done in lemma then the proof follows from 
computations that use property ( [M| , which itself uses all the couples (i,j) in 
V. However, the result of theorem!^ can not be applied directly to concrete net¬ 
works. This is the purpose of theorem which gives a sufficient condition for 
synchronization applicable in networks with arbitrary topology. Its proof relies 
on the computation of the sum of all minimal lengths joining any couple of nodes 
passing through a given edge of the graph. This idea comes from lemmaj^which 
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bounds all terms in the left hand side of (39) (in which all the couples (i,j) ap¬ 
pear) by terms of the right hand side (only the couples (i,j) corresponding to 
non zero coupling, i.e. Cij ^ 0, i.e. edges of the graph, appear). 

We start with the lemma that establishes the existence of a diagonal 
definite-positive matrix A that annihilates the effects of nonlinear terms. Let us 
denote by X an arbitrary vector of Xs the vector of the s first coordinates 
and Xd-s, the vector of the d — s last coordinates of X. 


Lemma 4. For each solution {Ui,U 2 , ■■■,Un) of (37), there exists a time T, a 
positive constant k and a definite positive diagonal matrix x M®, such 

that for all t>T, for all X S and for all i,j G {!,...,n}, 




[ DF{0Uj + {l-0)U^)de 
Jo 


X+ 


Xds- (^-a{x)Xd-s + [/o D<l>{0Uj + (1 - 0)u,)d0'\ < -K\\Xd- 

Proof. This follows from the L°°-bound of solutions of ([^, and by using young 
inequality. □ 


We will now present the main result that furnishes sufficient conditions for 
identical synchronization. To this aim, we introduce the following notations: 
for all i,j G {1, ...,n}, = Uj - u*, = Vj - Vi,i,j G {1, •■■,«}, Cij = 

and a = maxjg{i_,.._s} An. 

Theorem 2. If we assume that 


i<j 


i<j 


(39) 


then the system (37) synchronizes in the sense of definition 


Proof. We split G into the sum of two symmetric and antisymmetric matrices, 
E = ieik),i,k G {!,...,n} and L = {6ik)A,k G {l,...,n}: 


G = E + L, 


(40) 


where, 


dik 


2 {^ik C-ki) . 


(41) 


One can easily check that both matrices E and L have zero row sums. As 
Wij = Uj — Ui, Zij = Vj — ViAG G {1, we obtain, 


f aJijt — E{Uj.,Vj') Fi^Ui^Vi^ -)- QXwij I — 1 .jk CikVJik djkWjk dikWik)^ 

1 Zijt = -a{x)zij-G<^{x,Uj)-<^{x,Ui). 


( 42 ) 
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Besides, 


/» 1 7 ■ ^ 1 

F(uj,vj)-F{u,,v^)= —F(0UJ + il-O)U^)d^= / DF{dUj + (1 - 0)U,)d0 

Jo L^o 

where DF is the s x d Jacobian matrix of F, and, 

^{x,Uj)—^{x,Ui)= ( -jz^{0uj + {1 — 6)ui)d9 = [ D^{6uj + {1 — 9)ui)d9 

Jo L^o 

Hence, we can write. 


{U,-Ui), 




[ DF{9Uj + (1 - 9)Ui)d9 
Jo 


(Uj — Ui) + Ql\WijF 


F djkWjk 


fc=l 


Zijt — (T^X^Zij “t“ 


f D^{9uj + (1 — 9)ui)d9 
Jo 


Wii. 


where i,j S {1, Now, let us introduce the following function: 


(43) 


+ \Fj\in) ■ 

i=i j=i 


(44) 


In order to reach the synchronization, it is sufficient to find conditions ensuring 
that H is a Lyapunov function with negative orbital derivative. As the graph 
is connected, it would be natural to include in V, only terms Wij corresponding 
to non-zero coefficients eij. However, as we will see, including all terms wij in 
V will allows to vanish sums. 

Then, 


2 dt 


71 TL p 

= T.T. I {Wtj.Wijt + Zij.Zijt) 
i=l j=l 

71 n p r pi 

= E E ( / / DF{9U, + (1 - 9)U.)d9 

i=ij=i J^ Lio 

71 

F ^ ^, i.^jk’aJjk CikWik F SjkWjk dik‘OJik)) 


{U, - UJ) + QAw, 


k=l 


+ f Ztj.(—a(x)zij + f D^{9uj + {1 — 9)ui)d9 

Jo, Jo 


Wij)). 


Now by lemma and using green formula, there exists a definite positive 
diagonal matrix A such that, 
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1 dV ^ 

2 dt — 


n n / „ „ / n \ 

V] V] / AWij.W.j + / Wij. ^{CjkWjk - f-ikWik) 

^=l \k=i J 

+ J (^^{SjkWjk - S^kWik)^ • 

The last term in the above equation vanishes, indeed: 

n 

^ ^ - jk ^ik'^ik^ — ‘^^^i^j^k—l •^jk'^ik (bcCH-llSG 1 ^jk — 0) 

u — 1 

~ '^k^j SjkWij.Wik- 


i,i,k—l 


Obviously since Sjj = 0, 








Moreover, since Sjk = —^kj, we have, for all i G {1, ...,n}, 

n 

EE ^jk’dJij-Wik — 0. 

j = l k^j 

Therefore, asymmetric connectivity matrices with zero column-sums and 
row-sums can be treated as symmetric matrices. Now we deal with the other 
terms. We have: 

n n n n n 

EE Wij-'^iCjkWjk - CikWik) = y] CjkWij.Wjk- y] CikWij.Wik 

i—1 j — 1 k—1 i^j^k—l i^j^k—l 

n n 

— ^ ^ ^jk'^ij’’^jk ^ ^ ^jk'^ji-'^jk 

i,j,k—l i,j,k—l 

n 

— 2 ^ ^ CjkWij.Wjk 

ij,k—l 

n 

— 2 ^ ^ Cjk{Wik.Wjk ^kj^f 

i^j^k—1 


since Wij = Wik + Wkj- Moreover, 


n 

W^k-CjkWjk < Y 2^jk{\\Wik\\^ -G \\Wjk\\^) 


iJ^k—1 


ij,k—l 
1 
2 


j = l j,k=l 


i,k—l 


= ^ E ^JkWwjkW^, 


3,k=l 
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since ejk = 0. Finally, we obtain: 


i=i 


n „ „ n 

2lW — + / Awij.Wij) — n / ( ) 

i,i = l j.fe=l 

= 2V(-K|zij|^_j^+ / Aw,yWij)-2n I (y'(e„||w„|p) 

■In -In 


< 


z<j 


^ijl2,n + |w 

i<j 


ij l 2 ,a 


- 2 n^ 


Cij\Wij\2p 


i<3 


< 0 , 


if (wij, Zij) ^ 0 thanks to hypothesis 


□ 


Now we would like to prove theorem which gives a sufficient condition 
for synchronization which can be applied for general networks. An apparent 
difficulty comes from the fact that in the left-hand side of (39), all the variables 
appear while in the right-hand side, because of the parameters Cij , only variables 
corresponding to edges appear. The following lemma allows to obtain a bound 
with only terms corresponding to edges. 

Lemma 5. For all i,j G {l,...,n}, and for all sequence, (*z)/g{o,...,fc}: with 
io = i ,..., ik = j, we have: 


k-l 




(45) 


;=o 


Proof. We write: 


k-l 


1^0 
k-l 


k-l 


W. 


nn+i) 


k-l 


Z =0 






z=o 

k-l 


l<m 


< ^Ei 


m 


*i*! + l I 


by Young inequality. 




□ 

Note that the lemma is valid for an arbitrary sequence, but we will use it 
for sequences corresponding to edges in the graph. 
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Without taking care of the edges direction in the graph, for each {i,j),i<j, 
we choose a unique path of minimal length in the graph joining nodes i and j. 
We denote this path Pij, l{Pij) its length, and its nodes by: {ii)i^{o^...,iength(Pij)}i 
with io = i, ■■■,iiength{Pij) = j- For each (fc,l), k < I corresponding to an edge 
in the graph (i.e. Cki 7 ^ 0 ), we define a^i as the sum of all the lengths of the 
minimal paths passing trough the edge (fc, Z). If {k,l) is not an edge of the 
graph, we set otui = 0 : 


Oikl = 


J2i<j,{k,i)ePij if is an edge of the graph, 

0 if {k,l) is not an edge of the graph. 


(46) 


Figure gives an example of the computation of these coefficients. 
Theorem 3. Let us assume that for each edge (fc, 1), we have, 


-Oikl < CkU 
n 


then system (371 synchronizes in the sense of definition 

Proof. Now, the proof of theorem follows from theorem and lemma We 
have, 


El 

i<j 


W. 


length{Pij) — l 

< length{Pij) thanks to lemma 

i<j 1=0 

= ^ akiW'WkiW'^ by reordering the terms of the sum along the edges of the graph. 

k<l 

Now the result follows from theorem]^ Such a result has been found in 00 
for ODE systems. □ 



Figure 1: Example of calculus of the number aki- The minimal path 
joining 1 and 2 is P 12 = 1 — 2, P 13 = 1 — 2 — 3, P 14 = 1 — 2 — 4, 
Pi 5 = 1-2-4-5, P 23 = 2-3, P 24 = 2-4,P25 = 2-4-5, P 34 = 3-4, 
P 35 = 3 — 4 — 5, P 45 =4 — 5. It follows that 023 = 2 + 1 = 3. 


Corollary 1 (Fully connected network). We assume thatyi,j G n}, 

i ^ j Cj > f,, then (37) synchronizes in the sense of definition^ 


Proof. This result comes obviously from theorem since for fully connected 
network, aij = 1. Note that in the particular case of fully connected network, 
we could also conclude from theorem [2l □ 
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Corollary 2 (Unidirectionally ring network). We assume that the connectivity 
matrix G = {cij),l < i,j < n, is given hy cu = —c < 0 Vz G {l,...,n}, 
Cii+i = c„i = c\/i G {1, ...,n — 1}, and Cij = 0 otherwise. Then, if we assume 


that 


c > 


f^in^ + 8) 


*/ 

if 

if 


n odd 
n and - 
n even 


and ^ 


odd 


the network (37) synchronizes in the sense of definition^ 


Proof. We start with the case of n odd. For each couple of nodes of the graph, 
there is a unique path of minimal length joining the nodes. If we suppose that 
n = 2k 1, then for each node indexed by 1: 


Qfji+l — (1 ■■■ + fc) + (2 + — + /c) + ... + (fc — 1 + fc) + /c 

= V* i 

A^i—1 A^j—i J 

_ (n—l)(n+l)n 

“ 24 * 

The figure l^a gives an example of such a network. In the case where of n even, 
we assume that n = 2k. For each couple of nodes (z,j) in the graph, if the 
distance between i and j is less than fc, there exists a unique path of minimal 
length joining the nodes. But if the distance between i and j is equal to fc, there 
exists two distinct paths of minimal length joining the nodes. Therefore, we 
can choose for each couple {i,j) of distance k, alternatively the minimal path 
trough the left and trough the right (i.e. for example, for node 1 , the minimal 
path of length k trough the left). Then we find, if ^ is even. 


— (f T ■■■ “t“ ^) “t“ (2 + — k — 1) T ... + (/l — 1 k') k — 1 

— V'" 1 - kk 

~ 2 '^ 

_ ra(n^+2) 

~ 24 

Figure 0b gives an example of such a network. If ^ is odd and in the worst 
case, 

= (1 k) (2 k — 1) + (3 + — k) {k — 1) k 
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4 Numerical simulations 


In this section, we consider networks of type (37) with: 

d = 2, s = 1, 

F{u, v) = -{—u^ -\-3u — v),Q = — 
cr{x) = b, <f>(x, u) = au-\- c, 
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9 

P 


Figure 2: Unidirectionally ring. In panel a), the graph has an odd number 
of nodes, n = 7. There exists a unique minimal path joining each couple of 
nodes in the graph. For example, the computation of ai 2 is given by ai 2 = 
(1 + 2 + 3) + (2 + 3) + (3). The order of the computation follows from the 
counting of all the lengths of the minimal paths passing trough the edge ( 1 , 2 ), 
starting at 1,7 and 6 . In panel b), the graph has an even number of nodes, 
n = 8 . If the distance between two nodes is equal to ^ = 4, there exist two 
distinct minimal paths joining these nodes. For example, we can link node 1 
and node 5 either by the path 1 — 2 — 3 — 4 — 5 or by 1 — 8 — 7 — 6 — 5. Therefore, 
we choose the path 1 — 2 — 3 — 4 — 5 to connect nodes 1 and 5, whereas we choose 
the path 2 — 1 — 8 — 7 — 6 to connect 2 and 6 , and so on. Then computation of 
ai 2 is given by ai 2 = (1 + 2 + 3 + 4) + (2 + 3) + (3 + 4). 

i.e. each node is represented by the following reaction-diffusion system of 
FitzHugh-Nagumo type, 



(47) 


where u = u(x,t),v = v(x,t), du,a,b > 0, Q C is a regular bounded open 


set, and with Neumann zero flux conditions on the boundary. 

We perform numerical simulations in two cases, fully connected networks in 
one hand and unidirectionally coupled ring networks in the other hand, both 
with linear coupling. As far as we know, only few results exists for networks 
of reaction-diffusion systems, the case of a chain network has been partially 
investigated in [2]. We use the following parameters values: 


(48) 


0 = 1,6 = 0.001, e = 0.1, du = 0.05. 


The numerical integration of the network is realized using a C program, on 
n X [0, T] = [0,100] X [0,100] X [0, 3000]. The main concern of this section are: 

• heuristic laws of strength coupling with respect to the number of nodes in 
the graph 

• influence of space heterogeneity in initial conditions with regard to the 
synchronization phenomenon. 

Our main conclusion is that heterogeneity in initial conditions does not af¬ 
fect the general law of synchronization with respect to the number of nodes in 
the graph. As for ODE’s, and in good agreement with our theoretical results, 
the threshold value for synchronization is given by a ^ law in case of a fully 
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connected network, see figure]^ and a law in case of a ring network, see figure 


12 . However, it appears that the space heterogeneity of initial conditions will 
increase the threshold value, see figure Besides, the persistence of patterns 
require some symmetry. This kind of question have been mentioned for example 
in [3 EH, and the persistence of patterns, with respect to the size of fl and the 
diffusion coefficient have been studied in Hi- If there is not a certain kind of 
symmetry the network evolves towards homogeneous solutions whereas hetero¬ 
geneous patterns persist under some symmetry conditions on initial conditions. 
In the case of one node, this question is theoretically treated in a forthcoming 
paper [3] . Indeed, in figure we show how the increase of coupling strength 
for non symmetric initial conditions, make the patterns disappear. Most of our 
simulations are done with c(x) = 0 which correspond to the periodic regime of 
the FHN system. However, we also perform some numerical simulations with 
c(x) depending on x. This allows us to take advantage of both excitability and 
oscillatory regime of the FHN system. The emerging laws are the same in this 
case. 

We split this section in two subsections corresponding to the fully and uni¬ 
directional connected ring networks. 


4.1 Fully connected network 

For the fully connected network, the system reads as: 


euit = -Ui 3ui - Vi -f duAui - gn ^ {ui - Uj) 
Vit = aui - bvi + c{x) 


1 = 1, n, (49) 


where gn is the constant coupling strength between each couple of nodes in the 
graph. We present here five figures. 

• The figure represents the evolution of the network, for n = 3, when 
we choose two initial conditions with spirals and one homogeneous. The 
simulations show that we obtain synchronization for a threshold value of 
(73 ~ 0.015. They also show that we obtain asymptotically three spirals 
patterns. 

• In figure we represent the evolution of the coupling strength value gn 
needed for synchronization with respect to the number of nodes in the 
graph. The number of nodes is varied from 3 to 20, and we obtain an 
heuristic law in -. This is in good accordance with our theoretical results 
and is similar with previous results for ODE’s networks. We have also 
obtained the same heuristic laws when we choose only spirals for initial 
conditions or only uniform laws on [—1,1] or even with c{x) non-constant 
function. Thus the spatial effects induced by initial conditions do not 
affect this heuristic law already found on ODE’s fully linearly coupled 
networks. 
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• In figure we show a possible effect of increasing the coupling strength. 
For some values of the coupling strength, if the initial conditions are near 
a symmetric configuration a pattern will persist otherwise they asymp¬ 
totically evolve to a space homogeneous periodic solution, see and the 
references therein for a discussion in the case of a single reaction-diffusion 
system of FHN type. Here in the case of the fully connected network of 
FHN systems, for a coupling strength = 0.5, and with to spirals and one 
homogeneous for initial conditions, we observe that the network evolves 
toward a space homogeneous synchronized state: the system synchronizes 
but there is no more patterns. 

• The figure is obtained in an analogous way as the figure but with 
different initial conditions. Here, for each x G H, we choose Ui(x,0) and 
Vi(x, 0), i € {1,2,3} as random uniform values in [—1,1]. The simulations 
show that we obtain synchronization for a threshold value of gs ~ 0.014. 
They also show that we obtain asymptotically three similar patterns. 

• In figure as in figure we represent the evolution of the coupling 
strength value needed for synchronization with respect to the number 
of nodes in the graph. The number of nodes is varied from 3 to 20, and 
we obtain an heuristic law in ^. Thus the spatial effects induced by initial 
conditions do not affect this heuristic law. 


• In figure]^ we show the evolution of the coupling strength value g 2 o needed 
for synchronization with respect to the ratio between random uniform 
initial conditions and space homogeneous initial conditions in the fully 
connected network. The figure shows that increasing the space dispersion 
in initial conditions increase the coupling strength threshold. 


Finally, in we show the synchronization of a fully linearly connected 
network of type (491 with n = 3, and c(x) a non constant function: c{x) = 
0 if X is in a small neighborhood of (0,0), and c(x) = —I.l overwise. 
For one subsystem, and because of the oscillatory-excitable of the FHN 
system, this induces generation of periodic pulses starting at (0,0). For 
initial conditions, we choose three distinct space-homogeneous values. We 
observe that the synchronization occurs for gs > 0.013.. Asymptotically, 
the three subsystems evolve with spiral patterns. 
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Figure 3: Synchronization of a fully linearly connected network of type (491 
with n = 3, c(x) = 0. For initial conditions, we choose two spirals and one 
homogeneous in space. We observe that the synchronization occurs for > 
0.015.. Indeed, each row, from up to down, correspond respectively to the 
following values for g^: 0.005, 0.01, 0.013, 0.015. The two first columns represent 
the synchronization error respectively between ui and U 2 and U 2 and U 3 . The 
three others columns represent respectively from left to right the isovalues of 
ui, U 2 and U 3 for all x = {xi,X 2 ) € 11 at time 3000. Asymptotically, the three 
subsystems evolve with spiral patterns. 
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xlO"^ 



Figure 4: Evolution of the coupling strength value gn needed for 
synchronization with the respect to the number of nodes, n in the 
graph for a fully connected network of type (49) with c(x) = 0, when 
n is varied from 3 to 20. For initial conditions, we choose approx¬ 
imately 50% of spirals and 50% of homogeneous. The blue points 
represent the values obtained numerically, the red curve represents 


the function g„ = - 0.00043. Thus, we obtain heuristically a 

n 

i law which is common in other areas and highlights the synchro¬ 
nization emergent property, see for example [SlIlolIII]. 



0 500 10001500200025003000 0 500 10001500200025003000 0 20 40 60 80 100 0 20 40 60 80 100 

Time Time xl xl 


(a) (b) (c) (d) 


Figure 5: Synchronization of a fully linearly connected network of 
type (49) with n = 3, c(x) = 0 and g^ = 0.5. For initial conditions, 
we choose two spirals and one homogeneous in space. This shows 
an effect of the increasing of the coupling strength. Indeed, the 
simulation is the same as in figurej^with a coupling strength equal to 
0.5. We observe that the spiral patterns disappear, and the network 
evolve to space homogeneous solutions. This comes from a lack of 
symmetry in initial solutions 
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Figure 6 : Synchronization of a fully linearly connected network of 
type (49) with n = 3, c(x) = 0. For initial conditions, we choose ran¬ 
dom initial conditions as follows: for each x € fl, we whose Ui(x,0) 
and Vi(x,0), i G {1,2,3} as random uniform values in [—1,1], We 
observe that the synchronization occurs for > 0.014.. Indeed, 
each row, from up to down, correspond respectively to the following 
values of g^: 0.001,0.005,0.01,0.014. The two first columns repre¬ 
sent the synchronization error respectively between ui and U 2 and 
U 2 and U 3 . The three others columns represent respectively from 
left to right the isovalues of Ui, U 2 and M 3 for all x = {xi,X 2 ) G 11 
at time 3000. Asymptotically, the three subsystems evolve with the 
same patterns. 
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Figure 7: Evolution of the coupling strength value gn needed for 
synchronization with the respect to the number of nodes, n in the 
graph for a fully connected network of type (491 with c(x) = 0, 
when n is varied from 3 to 20. For initial conditions, we choose 
random uniform values in [—1,1]. The blue points represent the 

values obtained numerically, the red curve represents the function 
0.042 

= -. Thus, we obtain heuristically a - law. Note that this 

.n . . , " 

heuristically law is independent of the spatial structure of the initial 

conditions. 
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Percentage of the random initial conditons 
in the fully connected network of 20 neurons 


Figure 8: Evolution of the coupling strength value g 2 o needed 

for synchronization with the respect to the ratio between random 
uniform initial conditions and space homogeneous initial conditions 
in the fully connected network. The blue points represent the 
values obtained numerically, the red curve represents the function 
g^{p) = 352 X 10"V - 11 X 10"® + 0.00305. 
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Figure 9: Synchronization of a fully linearly connected network of 
type (491 with n = 3, and c(x) = 0 if a; is in a small neighborhood 
of the center of c(x) = —1.1 overwise. For one subsystem, and 
because of the oscillatory-excitable of the FHN systeclearpagem,this 
induces generation o periodic pulses starting at (0,0). For initial 
conditions, we choose three distinct space-homogeneous values. We 
observe that the synchronization occurs for gs > 0.013. Indeed, 
each row, from up to down, correspond respectively to the following 
values for g^: 0.01 and 0.013. The two first columns represent the 
synchronization error respectively between ui and U 2 and U 2 and 
U 3 . The three others columns represent respectively from left to 
right the isovalues of ui, U2 and M 3 for all x = {xi,X2) S 11 at 
time 3000. Asymptotically, the three subsystems evolve with spiral 
patterns. 


4.2 Unidirectionally coupled ring network 

For the unidirectionally ring connected network, the system reads as: 

{ n 

euit = duAui - Ui3ui - Vi - gn (Mi - m^+i) _ 

z = l,n, (50) 

Vit = aui - bvi c(x) 

where gn is the constant coupling strength between each couple of nodes corre¬ 
sponding to an edge in the graph. As for fully connected networks We present 
here three figures. 

• The figure represent the evolution of the network, for n = 3, when 
we choose all initial conditions with spirals. The simulations show that 
we obtain synchronization for a threshold value of 33 ~ 0.001. They also 
show that we obtain asymptotically three spirals patterns. 

• The figure is obtained in an analogous way as the figure [T0| but with 
different initial conditions. Here, for each x S O, we choose Ui(x,0) and 
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Vi(x, 0), * S {1,2,3} as random uniform values in [—1,1]. The simulations 
show that we obtain synchronization for a threshold value of gs ~ 0.02. 
They also show that we obtain asymptotically three similar patterns. 

• In figure |12[ we represent the evolution of the coupling strength value 
gn needed for synchronization with respect to the number of nodes in 
the graph. The number of nodes is varied from 3 to 20, and we obtain an 
heuristic law in . This is in good accordance with our theoretical results 
and with previous results for ODE’s networks. We have also obtained the 
same heuristic laws when we choose only spirals for initial conditions or 
spirals and homogeneous or even with c(x) non-constant function. Thus 
the spatial effects induced by initial conditions do not affect this heuristic 
law already found on ODE’s ring linearly coupled networks. 
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Figure 10: Synchronization of a fully linearly connected network 
of type (50) with n = 3, c{x) = 0. For initial conditions, we 
choose three spirals. We observe that the synchronization occurs 
foi" 93 ^ 0.001. Indeed, each row, from up to down, correspond re¬ 
spectively to the following values for g^: 0.0003,0.0005,0.001. The 
two first columns represent the synchronization error respectively 
between ui and U 2 and U 2 and u^. The three others columns rep¬ 
resent respectively from left to right the isovalues of wi, U 2 and U 3 
for all X = {xi^X 2 ) G Q at time 3000. Asymptotically, the three 
subsystems evolve with spiral patterns. 
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Figure 11: Synchronization of a fully linearly connected network of 
type (501 with n = 3, c(x) = 0. For initial conditions, we choose 
random initial conditions as follows: for each x € fl, we whose 
Ui(x, 0) and Vi(x, 0), i G {1, 2, 3} as random uniform values in [—1,1]. 
We observe that the synchronization occurs for gs > 0.02. Indeed, 
each row, from up to down, correspond respectively to the following 
values of 33 : 0.001,0.005,0.02. The two first columns represent the 
synchronization error respectively between ui and U 2 and U 2 and 
M 3 . The three others columns represent respectively from left to 
right the isovalues of mi, M 2 and M 3 for all x = (xi,X 2 ) € at time 
3000. Asymptotically, the three subsystems evolve with the same 
patterns. 
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Figure 12: Evolution of the coupling strength value gn needed for synchroniza¬ 
tion with the respect to the number of nodes, n in the graph of the unidirection- 
ally ring connected network of type (501 with c(x} = 0, when n is varied from 3 
to 20. For initial conditions, we choose approximately 50% of spirals and 50% 
of homogeneous. The blue points represent the values obtained numerically, the 
red curve represents the function = 0.0000167n^ -I- 0.00062n -I- 0.02. 


5 Conclusion 

In this paper, we have considered a network of n coupled reaction-diffusion sys¬ 
tems. We obtained three main contributions. First, we prove of the existence 
of the network attractor, and therefore, within this attractor, we have analyzed 
the synchronization behavior. We found out theoretically threshold values for 
synchronization for general class of networks with linear coupling. Finally, we 
have performed numerical simulations with different kind of initial conditions 
and studied numerically the spatial effects on synchronization and pattern for¬ 
mation. Thanks for the numerical experiments, we exhibited heuristic laws with 
regard to the number of nodes in the graph. The general form of these heuris¬ 
tic laws does not depend on spatial heterogeneity of asymptotic behavior. Our 
numerical simulations also show the persistence of asymptotic spatial patterns 
in complex networks. In future work, we aim to apply these results to more 
realistic neuronal networks. 

Appendix 

Theorem 4. Let g,h and y G L^(R) three positive functions. We suppose that 
for all t > to-' 

1 -™ + '* 

and, 

/ i+r 

g{s)ds < oi, / h{s)ds < a^, / y{s)ds < 03 , (52) 

where r,ai, 02 , are positive constants. Then 


y{t-\-r) < (aar + a 2 )e“SVt > to- 


(53) 
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Proof. Let sq > ^o- By (51), we have: 




d , - p 
dt 


(54) 


We integrate (54) between t and Sq + r for t S [sq, sq + ^]- We obtain: 

^so+r pt' 


s +r t pSQ-rT pt 

-/,“ + r) _ e ~-^=0 < / exp{— / g{s)ds}h{t')dt'. 

\J t «/ Sq 

(55) 


It follows that: 


^ pSQ-<rr 

e~°‘^y{so + r) — e~ -^‘o < J h{t')dt'. 

By multiplying by exp{ g{s)ds}^ we find: 

g/sQ ®(®)e““^y(so + r) — y{t) < e^-o h{t')dt'. 


Then we integrate between sq and sq + ?' with respect to t. This gives: 

e~°‘^y{so+r) <03 + a 2 r, 

thus: 

y{so + r) < {as + a 2 r)e“L 


(56) 

(57) 

(58) 

(59) 

□ 


Corollary 3. Let y,h G L^(IR) two positive functions. We assume that fort 
t>to: 

nil 

(60) 


^<h 

dt 


and, 


/ t+r i‘t-\-r 

h{s)ds < 02 , / ds < 03 


where r, 02,03 are positive constants. Then 

03 


y{t + r)< -1- 02 , Vt > to- 


(61) 


(62) 


Proof. It follows obviously from theorem with oi = 0. We can also give a 
direct proof: 

y{so + r) - y{t) < J^Sq + rhds. (63) 


Then we integrate between sq et Sq + r with respect to t. We obtain, 

y{so + r) < — + 02 . 

r 


(64) 

□ 




Complex network of reaction-diffusion equations 


32 


References 

[1] B. Ambrosio and J-P. FRANgoiSE, Propagation of Bursting Oscillations, 
Phil. Trans. R. Soc. A, 367 (2009), pp. 4863-4875. 

[2] B. Ambrosio and M.A. Aziz-Alaoui, Synchronization and control of 
coupled reaction-diffusion systems of the FitzHuqh-Naqumo-type, Comput. 
Math. AppL, 64 (2012), pp. 934-943. 

[3] B. Ambrosio, M.A. Aziz-Alaoui, On the basin of attraction of patterns 
for the reaction-diffusion system of generalized FitzHugh-Nagumo type in 
ocillatory regime, in progress. 

[4] B. Ambrosio, Propagation d’ondes dans un milieu excitable: simulations 
numeriques et approche analytique, Thesis of University Pierre et Marie 
Curie-Paris 6, 2009. 

[5] M.A. Aziz-Alaoui, Synchronization of Chaos, Encyclopedia of Mathemat¬ 
ical Physics, Elsevier, Vol. 5, 2006 pp. 213-226. 

[6] V.N. Belykh, I. Belykh, M. Hasler, Connection graph stability method 
for synchronized coupled chaotic systems, Physica D, 195 (2004), pp. 159- 
187. 

[7] I. Belykh, V.N Belykh, M. Hasler, Sychronization in asymmetrically 
coupled networks with node balance. Chaos 16, (2006), 015102. 

[8] E. Conway, D. Hoff, J. Smoller, Large time behavior of solutions of 
systems of nonlinear reaction-diffusion equations, SIAM J. Appl. Math. 35 
(1978), pp. 1-16. 

[9] N. Corson, M.A. Aziz-Alaoui, Asymptotic dynamics of the slow-fast 
Hindmarsh-Rose neuronal system, Dyn. Contin. Discrete Impuls. Syst. Ser. 
B, 16 (2009), pp: 535-549. 

[10] N. Corson, M.A. Aziz-Alaoui, R. Ghnemat, S. Balev, C. 
Bertelle, Modeling the Dynamics of Complex Interaction Systems: from 
Morphogenesis to Control, Internal. J. Bifur. Chaos, 22 (2012), pp: 1250025. 

[11] N. Corson, M.A. Aziz-Alaoui, Complex emergent properties in synchro¬ 
nized neuronal oscillations, in M.A. Aziz-Alaoui and C. Bertelle (eds.): From 
System Complexity to Emergent Properties, Springer, (2009), pp. 243-259. 

[12] P. Dorfler and E. Bullo, On the Critical Coupling for Kuramoto Os¬ 
cillators, SIAM J. Appl. Dyn. Syst., 10 (2011), pp. 10701099. 

[13] P. Dorfler and P. Bullo, Synchronization and Transient Stability in 
Power Networks and Nonuniform Kuramoto Oscillators SIAM J. Control 
Optim., 50 (2009), p. 16161642. 



Complex network of reaction-diffusion equations 


33 


[14] R. A. FitzHugh, Impulses and physiological states in theoretical models 
of nerve membrane, Biophysical Journal. 1 (1961), pp. 445-466. 

[15] A. Friedman, Partial Differential Equations of Parabolic Type, Dover Edi¬ 
tion, 2008. 

[16] M. Golubitsky and I. Stewart, The Symmetry Perspective 
,Birkhauser, 2002 

[17] D. Henry, Geometric Theory of Semilinear Parabolic Equations, Springer, 
1981. 

[18] A.L. Hodgkin and A.F. Huxley, A quantitative description of mem¬ 
brane current and its application to conduction and excitation in nerve, J. 
Physiol.117, (1952), pp. 500-544. 

[19] O.A. Ladyzenskaja, V.A. Solonnikov and N.N. Ural’ceva, Linear 
and Quasilinear Equations of Parabolic Type, Am. Math. Soc., Providence, 
Rhode Island, Transl. of Math. Monographs 23, 1968. 

[20] J. L. Lions, Quelques methodes de resolution des problemes aux limites 
non lineaires, Dunod, Paris, 1969. 

[21] M. Marion, Einite Dimensionnal Attractors associated with Partly Dis¬ 
sipative Reaction-Diffusion Systems, SIAM, J. Math. Anal. 20 (1989), pp. 
816-844. 

[22] X. Mora, Semilinear Parabolic Problems define semiflows on spaces. 
Transactions of the AMS, 278 ( 1983), pp. 21-54. 

[23] J. Nagumo, S. Arimoto and S. Yoshizawa, An active pulse transmis¬ 
sion line simulating nerve axon, Proc. IRE. 50 (1962), pp. 2061-2070. 

[24] A. PiKOVSKY, M. Rosenblum and j. Kuths, Synchronization, Cam¬ 
bridge University Press, 2003. 

[25] J. Robinson, Infinite-Dimensional Systems, Cambridge University Press, 
2001 . 

[26] J. Rauch and J. Smolder, Qualitative Theory of the PitzHugh Nagumo 
Equations, Advances in Mathematics, 27 (1978), pp. 12-44. 

[27] E. Rothe, Global Solutions to Reaction-Diffusion Systems, Springer- 
Verlag, Berlin, 1984. 

[28] J. Smolder, Shock Waves and Reaction-Diffusion Equations, Springer, 
1994. 

[29] R. Temam, Infinite-Dimensional Dynamical Systems in Mechanics and 
Physics, Springer, 1988. 



